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1.  INTRODUCTION 

This  note  outlines  methods  which  are  described  in  detail  in  reference  1,  for 
calculating  the  curvature  characteristics  of  a reflected  electromagnetic  wave, 
given  the  curvature  characteristics  of  the  incident  wave  and  a geometrical 
description  of  the  reflecting  surface  at  the  point  of  reflection.  The  methods 
are  based  on  geometrical  optics  and  are  of  particular  use  in  estimating  double- 
bounce radar  cross  sections,  for  which  individual  radii  of  curvature  of  reflected 
waves  are  required,  when  the  radii  of  curvature  of  the  incident  wave  for  the 
second  reflecting  surface  are  finite.  This  work  has  application  in  the  calcula- 
tion of  the  radar  cross  section  of  an  aircraft. 

Since  the  standard  textbooks  on  radar  cross  section  estimation  were  published, 
simpler  methods  than  Fock's  tensor  method, (ref .3,6)  have  been  presented  for 
determining  the  curvature  characteristics  of  reflected  waves  for  the  general 
case, (ref. 1,7,8) . The  prime  reference  for  the  simpler  method  is  that  of 
Deschamps, (ref . 1)  with  application  of  the  formulae  for  specific  situations  being 
outlined  by  Kouyoumjian  and  Pathak(ref .7) , and  Lee, (ref. 8). 

Although  simple  examples  of  Deschamps  general  "mirror"  relations,  (which  are 
the  basis  of  the  method  outlined  below),  appear  in  other  references,  either 
explicitly, (ref .9)  or  implicitly, (ref. 5) , it  was  felt  that  their  usefulness  and 
applicability  were  not  widely  enough  appreciated.  This  document  aims  to  widen 
this  appreciation. 

Use  of  the  "mirror"  relations  of  Deschamps  requires  knowledge  of  the  principal 
radii  of  curvature  of  the  reflecting  surface  at  the  reflection  point  and  of  the 
orientation  of  the  principal  directions.  In  specific  instances,  such  as 
surfaces  of  revolution,  these  may  be  readily  determined  sometimes  by  inspection 
(ref. 2).  However  for  a general  point  on  a general  body  such  as  an  ellipsoid 
(which  is  frequently  employed  for  the  local  representation  of  aircraft  surfaces) , 
use  of  relations  from  three  dimensional  co-ordinate  geometry  is  required(ref . 10) . 
Appropriate  formulae  which  are  set  down  below,  have  been  programmed  on  a digital 
computer (ref . 11) , for  ease  of  solution  of  the  magnitude  of  the  principal  radii  of 
curvature  and  principal  directions  for  a specified  point  on  a given  ellipsoid. 

The  methods  outlined  below  have  been  verified  in  part(ref.2)  by  applying  them 
to  a number  of  double-bounce  situations  considered  in  the  textbooks(ref .3) , (viz. 
two  adjacent  spheres,  and  two  adjacent  paraboloids). 


2.  OUTLINE  OF  METHODS 

2.1  Radii  of  curvature  of  the  reflected  wave  for  simple  shapes 

It  is  convenient  to  describe  the  curvature  characteristics  of  a wave- 
front  by  its  curvature  matrix  Qr,  which  describes  the  way  in  which  the 
total  wavefront  curvature  is  divided  with  respect  to  a specified  axial 
system.  For  example,  when  the  two  principal  directions  of  curvature  of 
a wavefront  are  parallel  to  two  lateral  axes  of  the  specified  axial 
system,  the  curvature  matrix  is  diagonal,  with  the  diagonal  elements 
being  the  inverse  of  the  principal  radii  of  curvature  of  the  wavefront. 

In  a similar  way,  the  curvature  matrix,  CQ  for  a reflecting  surface, 

for  a given  reflection  point,  is  also  diagonal  when  referred  to  its 
principal  directions,  with  the  diagonal  elements  being  the  inverse  of  the 
local  principal  radii  of  curvature  of  the  surface. 

Clearly,  the  way  in  which  the  components  of  curvature  of  an  incident 
wave  are  modified  on  reflection  by  the  curvature  properties  of  the 
reflecting  surface  must  depend  on  the  relative  orientation  of  the 
respective  principal  directions  of  the  incident  wavefront  and  of  the 
reflecting  surface. 
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Deschamps  showed  that  the  curvature  matrix  Qr  of  a reflected  wave  is 

related (ref . 1)  to  the  curvature  matrix  of  the  incident  wave  and  to  the 

curvature  matrix  CQ  of  the  local  reflecting  surface,  to  the  angle  of 

incidence  0.  and  to  the  axial  transformation  matrix  G as  follows: 
i 

-i  -\T 


Qr  = C^1  ♦ 2 (G_I  )T  CQ  G'1  cos  0i 


where 


ni  A ri/p,i  0 -1 

^ L o i/p^J* 

Pi1,  Pi  1 being  the  principal  radii  of  curvature  of  the  incident  wave, 

ri/Ri  0 ~| 


co  = L o 


l/R.  J 


Ri , Rj  being  the  principal  radii  of  curvature  of  the  local  reflecting 
surface,  and  where 


* |V  * a 

Lx? 1 . U. 


Xi1 


Xi 1,  )C3  1 being  the  principal  directions  of  the  incident  wave  and  Uj  , Ifc 
being  the  principal  directions  of  the  local  reflecting  surface.  The 
principal  radii  of  curvature  (Pir,  PjT)  of  the  reflected  wave  are  the 
reciprocals  of  the  principal  values  of  Q . 

2.1.1  Application  when  the  matrices  are  diagonal 

In  the  case  of  a wave  incident  on  a surface  with  incidence 
angle  0 ^ and  with  one  of  the  surface's  local  principal  directions 

lying  in  the  plane  of  incidence,  figure  1,  as  in  two  double-bounce 
examples  of  reference  3,  then  matrix  G becomes  as  follows: 


- P 0 1 

|_0  cos 


If  the  local  principal  radii  of  curvature  of  the  surface  are 
Ri , R2  , then 


“ 1/Pi 
0 


1/Pj 


[o 


0 -lri/R, 

l/cos  0^_  _ 0 


0 If  0 1 

1/B.JLo  1/C05  »4J 


- 3 - 


WSRL-0030-TR 


0 1/Pj 


2 cos  0 . 


where 


r a rl/p,r  0 1 

f = 

L 0 l/p,rJ 


Rj  cos  0 . 


when  Qr  is  diagonal 


i r , . i 2 cos  0 • 
1/Pi  = 1/Pi  + 1_ 


1*>>r  ■ 

l 


2.1.2  Application  when  the  matrices  are  non -diagonal 

In  the  more  general  case  when  the  various  principal  directions 
are  not  either  coplanar  with  or  normal  to  one  another,  the 
principal  radii  of  the  reflected  wave  are  derived  from  the 
reciprocals  of  the  eigen  values  of  Q1.  Specific  algebraic 

r p 

formulae  for  Pi  , and  p2  • and  for  the  associated  principal 
directions  are  given  by  Kouyoumjian  and  Pathakfref .7) . 

2.2  Radii  of  curvature  of  the  local  surface  of  an  ellipsoid 

Let  the  ellipsoid  be  described  by  the  equation 

x2  y2  z2 


and  the  point  on  the  surface  where  the  radii  of  curvature  are  required 
be  (xs,  ys,  zs). 

The  following  relations  are  taken  from  Bell(ref.lO) . If  the  length 
of  the  perpendicular  from  the  origin  to  the  tangent  plane  at  (x  , y , z ) 
is  p,  where  s s s 

/ x 2 + y 1 t z 2\"!s 
p i*-  c+)  ' 

then  the  direction  cosines  of  the  normal  to  the  ellipsoid  at  (x  , y , z ) 

s s s 


l = pxs/a2 
» ■ pys/b* 
n = ptg/c2 
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The  square  of  the  lengths  (ri ,r2 ) of  the  semi-axes  of  section  through 
the  origin  parallel  to  the  tangent  plane  at  (xg,  yg,  zg)  are  found  from 

the  roots  ri2.,  rj2  of  the  following  equationfref  .10) : 

WJl  jL\ 

Vl>V  a’c’  ■ a’bV 

''(■5  *7 


-r2 


L 


l2  + m2  + n2  ‘ ■ 0 


The  principal  radii  of  curvature  R|  ,,  Rj  of  the  surface  of  the  ellipsoid 
are (ref .10) : 

Ri  = r,2,/p  and  Ra  = r2  Vp 

The  direction  cosines  describing  the  two  principal  directions  (Xi  ,#ii  ,v2  ) 
and  (X2,p2,vj)  are  found  by  substituting  ri2  and  r22  in  turn  in  the  follow- 
ing relations(ref . 10) : 


M?-‘)  •:<?-) 


where 


Xf2  + /i±2  + Vj2  = 1 , i = 1,2. 


2.3  Calculation  of  reflected  electric  fields 

Once  the  two  principal  radii  of  curvature  of  a reflected  wave  (Pir,  P2 T) 

are  known,  the  relative  value  E of  the  electric  field  voltage  at  a 

r 

distance  R from  the  reflection  point  in  free  space  for  unity  value  of  the 
indicent  electric  field  voltage  is  obtained (ref .3,4,5)  from: 


E 

r 

E. 

1 


r p»r'  Pi*  1 h 

l(Pir.  + R)  (p2r  ♦ R)J 


with  E.  = 1. 

1 


2.3.1  Use  in  RCS  calculation 

The  above  relation  for  E / E.  can  be  used  to  obtain  the 

I r x 1 

resultant  radar  cross  section  (0)  as  R ♦ 00  from 


a 


Limlt  4»B?  ' 

R ♦ — 
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2.3.2  Double  bounce  situations 

Alternatively,  if  double  reflections  are  involved,  the 
above  relation  for  |Er/Ej|  enables  determination  of  the 

magnitude  of  the  field  voltage  incident  at  the  second 
reflection  point  by  substituting  the  appropriate  range  R in 
the  above  expression. 

For  use  of  Deschamps'  relations  for  the  curvature  matrix  at 
the  second  reflection  point,  it  is  clearly  necessary  to  increase 
the  calculated  radii  of  curvature  of  the  first  reflected  ray  by 
the  separation  of  the  two  reflection  points  to  obtain  the  radii 
of  curvature  of  the  second  incident  ray  at  the  second  reflection 
point, (ref .2) . 


3.  DISCUSSION 

In  most  scattering  problems  utilizing  geometric  optics,  a basic  computation 
step  is  the  determination  of  the  reflected  field  for  a specified  point  of 
reflection  on  a given  conducting  surface,  knowing  the  characteristics  of  the 
incident  wavefront. 

In  the  general  case,  (i.e.  when  the  principal  directions  of  the  scattering 
surface  are  not  necessarily  either  coplanar  with  or  normal  to  the  plane  of 
incidence),  the  standard  textbooks  on  radar  scattering  from  simple  shapes 
(ref. 3, 4, 5),  fall  back  on  Fock's  method(ref .6)  which  uses  tensor  notation  and 
was  originally  published  in  1950.  Reference  3 gives  a detailed  description 
and  examples  of  the  method.  However,  the  required  calculation  in  the  general 
case  is  more  readily  handled  by  Deschamps'  relations(ref .1) , using  the  more 
familiar  matrix  notation,  as  it  is  indeed  in  some  simple  situations  now  to  be 
discussed. 

As  noted  in  para  2.1.1  above,  when  one  of  the  surface's  principal  directions 
is  coplanar  with  the  plane  of  incidence  which  also  contains  a principal 
direction  of  the  incident  wavefront  the  relevant  matrix  becomes  diagonal  and 
the  curvatures  of  the  reflected  waves  are  given  by: 


and 


1/Pi 1 + 


2 cos  0 . 
k 

Ri 


1/Pa  r- 


1/P, 1 


2 

R,  cos  0 . 

i 


From  para  2.3  above,  at  a distance  R from  the  point  of  reflection,  the 
relative  magnitude  of  the  reflected  field  voltage  to  the  incident  field  voltage 
at  the  reflection  point  is  equal  to: 


[ 


r r 
Pi  Pa 

(Pi'  + R)  (Pa1  + R)J 


% 


{ (1  + R/Pi r.)  (1  ♦ R/Pa  *•) } 
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V 


r r 

Substitution  of  the  above  relations  for  1/Pi  , 1/Pj  gives: 


1 + — . + 2R  cos  fli\  /i  ♦ JL,  + — 

Pi  X'  Ri  • / \ Pj * Rj  cos  0 . 


In  the  case  of  a point  source  located  at  distance  d from  the  reflecting  point. 

Pi 1 = Pj 1 = d and  the  aoove  expression  then  becomes  identical  in  form,  (but  with 
different  notation),  to  that  quoted  by  Bowman  et  al(ref.5,  pg  24).  This 
latter  relation  is  used  extensively,  but  without  derivation  in  reference  5, 
since  many  of  the  scattering  examples  considered  there  relate  to  surfaces  of 
revolution  with  the  plane  of  incidence  containing  the  axis  of  revolution  and 
hence  containing  a principal  direction  of  the  surface  at  the  reflection  point. 

It  is  seen  that  Deschamps'  relations  readily  provide  the  necessary  relations  for 
such  simple  situations,  as  in  fact  they  do  for  more  complex  situations. 

The  geometrical  description  of  common  scattering  surfaces  also  differs  some- 
what between  authors  when  utilizing  methods  of  geometrical  optics  for  the 
solution  of  scattering  problems.  Thus  Bowman  et  al(ref.5)  utilize  oblate 
spheroidal  co-ordinates  rather  than  Cartesian  co-ordinates  when  describing  high 
frequency  scattering  from  a conducting  oblate  spheroid,  (an  ellipsoid  of 
revolution) . The  choice  of  such  co-ordinates  enables  the  high  frequency 
solution  to  be  considered  as  a special  case  of  the  solution  for  an  arbitrary 
frequency.  However,  when  prime  interest  is  in  a geometrical  optics  solution, 
the  use  of  Cartesian  co-ordinates  and  the  standard  methods  of  co-ordinate 
geometry(ref . 10)  would  appear  to  provide  the  simplest  methods  for  most  people. 

Although  the  techniques  outlined  above,  based  on  Deschamps'  relations,  are 
satisfactory  for  smooth  conductors  of  finite  (but  non-zero)  radii  of  curvature, 
there  is  also  a need  in  practice  for  similar  simple  techniques  by  which  to 
describe  specular  reflections  from  bodies  such  as  finite  cylinders  of  length  L, 
having  one  infinite  radius  of  curvature. 


4.  CONCLUSIONS 

Deschamps (ref . 1)  has  provided  a simple  matrix  relation  by  which  to  compute 
the  components  of  curvature  of  short  electromagnetic  waves  following  their 
reflection  from  conductors  of  simple  shape. 

Such  information  directly  enables  calculation  of  the  relative  magnitude  of 
the  specularly  reflected  field  at  a given  distance  from  a simple  conductor, 
provided  its  principal  radii  of  curvature  at  the  reflection  point  are  finite 
but  non-zero.  The  relation  used  twice  in  sequence  readily  enables  the 
calculation  of  double-bounce  radar  cross  sections  for  two  adjacent  simple  shape 
conductors . 

Deschamps'  relations  were  published  subsequent  to  the  standard  books  on  radar 
cross  section  estimation  (ref  .3,4,5)  which  qu^.-.e  the  more  complicated  tensor 
relations  published  by  Fock  in  1950  to  meet  the  same  requirement (ref . 6) . 

Because  of  the  generally  greater  familiarity  with  matrices  than  tensors, 
Deschamps'  relations  present  a significant  advance  in  means  of  estimating  the 
characteristics  of  reflected  waves. 
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NOTATION 

Curvature  matrix  of  a local  reflecting  surface 
Voltage  of  incident  field 
Voltage  of  reflected  field 

Transformation  matrix  connecting  the  principal  directions 
of  the  reflecting  surface  and  those  of  the  incident  wavefront 

Length  of  cylinder 

Matrix  describing  the  components  of  curvature  of  an  incident 
wavefront 

Matrix  describing  the  components  of  curvature  of  a reflected 
wavefront 

Range  from  reflecting  point 

Radii  of  curvature  of  a local  reflecting  surface 

Unit  vectors  describing  the  principal  directions  of  the  local 
reflecting  surface 

Unit  vectors  describing  the  principal  directions  of  the 
incident  wavefront 


Lengths  of  semi-axes  of  an  ellipsoid 

Separation  of  two  reflection  points,  or  distance  of  a point 
source  from  a reflection  point 

Direction  cosines  of  local  normal  to  an  ellipsoid 

Length  of  perpendicular  from  the  origin  to  a tangent  plane 
Lengths  of  semi  axes  of  an  elliptical  section 

General  Cartesian  co-ordinates 


Co-ordinates  of  a point  of  reflection 


\ 
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6. 


i 
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Angle  of  incidence 

Wavelength  of  electromagnetic  radiation 

Direction  cosines  of  the  principal  directions  of  a 
local  reflecting  surface,  i = 1,  2 

Principal  radii  of  curvature  of  an  incident  wavefront 


Principal  radii  of  curvature  of  a reflected  wavefront 


a 


Radar  cross  section 
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